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In this work, natural neutron spectra at nine sites in Tibet region were measured using a multi-sphere neutron 
spectrometer. The altitude-dependence of the spectra total fluence rate and ambient dose equivalent rate were 
analyzed. From the normalized natural neutron spectra at different altitudes, the spectrum fractions for neutrons 
of greater than 0.1 MeV do not differ obviously, while those of the thermal neutrons differ greatly from each 
other. The total fluence rate, effective dose rate and the ambient dose equivalent rate varied with the altitude 


according to an exponential law. 
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I. INTRODUCTION 


The issue of natural radiation effects on human health has 
always been a focus of the United Nations Scientific Com- 
mittee on the Effects of Atomic Radiation, and an important 
part of the issue is about the health effects of natural neu- 
trons [1]. Natural neutrons near ground can be produced from 
several modes, including the interaction of cosmic ray parti- 
cles with the atmosphere and surface medias, the capture of 
cosmic ray particles on the earth and the ground natural ra- 
diation [2], and the energy spectrum spans over ten orders of 
magnitude, from thermal neutrons to high energy neutrons of 
hundreds of MeV [1, 3]. The measurement of natural neutron 
energy spectrum is an international hotspot, and laboratories 
in USA, Japan and other countries have been doing so for half 
a century [4-11]. Especially, Japan has completed a nation- 
wide survey of the natural neutron radiation level [4]. Al- 
though many works on natural neutron dose have been done 
in China, there is still a need of measuring the nationwide 
spectrum of natural neutrons [10, 11]. Due to the high alti- 
tude and high cosmic ray radiation level in Tibet, the natural 
radiation dose is higher in Tibet than those of low altitude ar- 
eas. In this work, measurements with a multi-sphere neutron 
spectrometer were performed on natural neutron spectrum in 
areas of different altitudes in Tibet. Altitude variations of the 
spectrum, total fluence rate, and ambient dose equivalent rate 
were analyzed. The results provide technical supports and 
basic data for the radiation survey and radiation hazard as- 
sessment on natural neutrons in Tibet. 


Il. EXPERIMENTAL 
A. Instrument 
The multi-sphere neutron spectrometer used in this work 


is shown schematically in Fig. 1. A ®5.08cm *He spheri- 
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cal proportional counter (LND-27036) was used as the ther- 
mal neutron detector. The moderator was made of high- 
density polyethylene (0.94 g/cm?), being 8cm, 11cm, 15cm 
and 23 cm in diameter or nothing (bare detector). The moder- 
ator balls and the configuration are shown in Fig. 2. 


Polyethylene moderator 


3He proportional 
counter 


Fig. 1. (Color online) Schematics of the multi-sphere neutron spec- 
trometer system. 


Fig. 2. (Color online) The moderator balls and detector 
configuration. 
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Fig. 3. (Color online) The response functions of the multi-sphere 
neutron spectrometer. 
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Fig. 4. (Color online) Spectra of natural neutrons at spots of different altitudes in Tibet. (a) as-measured, (b) normalized. 


TABLE 1. The altitudes, longitudes and latitudes data 


Spot No. Altitude (m) Longitude Latitude Places 

1 1 946 95°0'52.1"E 30°1°45.2"N Lulang Town 

2 2677 95° 46'2.4”"E 29°51'26.8"N Bome County 

3 2955 94°21'39.9"E 29°38'1.9"N Bayi Town 

4 3 200 93°35'36.5”"E 29°8'36.4"N Gyaca County 

5 3 490 90°53'36.3”E 29°17'20.5"N Shannan Prefecture 

6 3 826 88°51'30.3”E 29°17'21.4"N Shigatse Prefecture 

7 4 000 89°5'33.5”E 29°40'23.7”"N Namling County 

8 4 544 90° 40'39.2”E 30° 16'2.6"N Nyenchen Tanglsha Mountains pass 


The response functions of the multi-sphere neutron spec- 
trometer were calculated with the MCNPX code, using the 
method in Ref. [12], and the results are shown in Fig. 3. Be- 
fore the survey, the neutron spectrometer was tested, so as 
to check its performances. Energy resolutions of the detec- 
tor coated with different moderator balls were measured at 
about 7% at 765 keV of reaction *He(n, p)T, while the detec- 
tion efficiencies were about 10% as calibrated with 7°?Cf and 
241 Am-Be neutron sources. These agreed well with the simu- 
lation results. 


B. Measurement 


Spots of different altitudes in Tibet were chosen and natural 
neutron spectra were measured at each of the spots using the 
neutron detector coated with different moderator balls. Con- 
ditions of the spectrometer system kept unchanged to measure 
all the neutron spectra. And each spectrum was recorded in 
over two hours, with the counts of each moderator ball being 
greater than 300. The Global Positioning System was used 
for getting the coordinates and altitudes of the measurement 
spot. 


Ill. RESULTS AND ANALYSES 
A. Neutron energy spectra 


The natural neutron spectra obtained using an unfolding 
program of the least square method are shown in Fig. 4(a), 


with the altitudes, longitudes and latitudes listed in Table 1. 
It can be seen that the neutron fluence rates of each en- 
ergy range increase with the altitude, because the cosmic- 
ray intensity increase with the altitude. Figure 4(b) shows 
the neutron spectra normalized to the total fluence rate. The 
spectrum fractions of neutrons of greater than 0.1 MeV do 
not differ obviously, while those of the thermal neutrons at 
different altitudes differ greatly from each other, due possibly 
to differences in ground media and vegetation at the measure- 
ment sites of different altitudes, hence the different modera- 
tion and absorption to the natural neutron. 


B. Neutron fluence rate and the estimated neutron dose rate 


Total fluence rate, effective dose rate and ambient dose 
equivalent rate were calculated from the measured spectra at 
each altitude, as given in Table 2. In the calculation, the ef- 
fective dose per unit fluence for monoenergetic neutrons in- 
cident in isotropic geometries to adult male and female ref- 
erence phantoms calculated by PHITS [13] were used, and 
the conversion coefficients of ambient dose equivalent rate 
for monoenergetic neutrons were adopted from the GBZ/T 
202-2007 [14]. 

By fitting the data in Table 2, the altitude variations of neu- 
tron fluence rate and neutron dose rate can be obtained, as 
shown in Fig. 5. It can be seen that the total neutron fluence 
rate, effective dose and the ambient dose equivalent rate in- 
crease exponentially with altitude of the measurement sites, 
as pointed out in Refs. [1, 15]. 
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TABLE 2. Calculated neutron fluence rate and dose rate at various altitudes in Tibet 


Altitude Atmospheric depth Fluence rate Effective dose rate Ambient dose equi- 
(m) (g/cm?) (cm7! 871) (ISO geometry) (nSv/h) valent Rate(nSv/h) 
1 946 811 0.0239 28.2 40.9 
2677 741 0.0388 48.4 70.1 
2955 716 0.0513 63.3 91.6 
3 200 695 0.0560 68.6 99.4 
3 490 671 0.0620 72.6 105.3 
3 826 644 0.0679 87.4 126.4 
4 000 630 0.0984 129.2 186.6 
4544 589 0.1254 156.5 226.6 
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Fig. 5. (Color online) Altitude dependence of the total neutron fluence rate (a), effective dose rate (b) and ambient dose equivalent rate (c). 


TABLE 3. Standard deviations (SD) of the a-values 


where y is in cm~? s71, his in m, and E and H * (10) are in 


nSv/h. 


Standard deviations of the a -values are given in Table 3. 


Parameters a (m+) x1074 SD (m7!) x 1074 
p 6.25 0.54 
E 6.42 0.64 
H x (10) 6.41 0.62 


TABLE 4. The a-values at various geomagnetic latitudes obtained 
by different groups 


Geomagnetic latitude (°) a-value (x107? m7") 


18-21 0.625*(This work) 
26 0.7 [7] 

48 0.78 [1] 

49 0.85 [16] 

50 1.04 [17] 


The variations of the total neutron fluence rate y, effective 
dose rate E and ambient dose equivalent rate H » (10) can be 
fitted by exponential function of R- e~, where, h is altitude 
of the measurement sites, R and a are the fitting parameters. 
The fitted functions were obtained as Eqs. (1), (2) and (3). 


y= 0.0073e°: 0006252 (1) 
E= 8.61e0-000642h (2) 
H * (10) = jeer (3) 


It can be seen that the a-values of y, E and H x (10) are 
almost the same. This means that the natural neutron spectra 
changed just a little, so all the a-values calculated from the 
spectra varied with the altitude nearly in the same way. 

Table 4 compared the a-values obtained at various geo- 
magnetic latitudes by different groups. Both the a -value and 
the latitudes (18°—21°) of the present work are the smallest. 
This confirms that the a-values increase with the geomag- 
netic latitude, and that magnitude of the attenuation of nat- 
ural neutrons is influenced by the geomagnetic intensity [1]. 
As aresult, in order to assess the natural neutron dose in Tibet 
more accurately, the local a-value should be applied. 


IV. CONCLUSION 


Natural neutron spectra were measured at nine spots of 
different altitude: (1946—-4544m) in Tibet region. The 
effective dose and ambient dose equivalent rate were calcu- 
lated. The following conclusions can be obtained: 

In the normalized natural neutron spectra at different al- 
titudes, the spectrum fractions of neutrons of greater than 
0.1 MeV do not differ obviously, but those of the thermal neu- 
trons differ greatly from each other; 

The natural neutron fluence rate of each energy range in- 
creases with the altitude by y = 0.0073e° 00854 cm~? s—!; 

The effective dose rate of natural neutrons varies with the 
altitude by Æ = 8.61e° 0006428 nSv/h. 
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